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Q\ \ The explosion of a white dwarf of mass 1.36 Mq has been simulated in three dimensions 

^ ■ with the aid of a SPH code. The explosion follows the delayed detonation paradigma. In 

O ' this case the deflagration-detonation transition is induced by the large corrugation of the 

^ . flame front resulting from Rayleigh- Taylor instability and turbulence. The nucleosynthetic 

Q ! yields have been calculated, showing that some neutronized isotopes such as ^^Fe or 

"Th I ^^Ni are not overproduced with respect to the solar system ratios. The distribution 

Q-i' of intermediate-mass elements is also compatible with the spectra of normal SNIa. The 

Q ■ exception is, however, the abundance of carbon and oxygen, which are overproduced. 

^: 

■ 1. Introduction 



The undeniable influence of Type la Supernovae explosions in, among others, the chem- 
ical evolution of the galaxy and cosmology makes the quest of solid theoretical models 
5^ I as an urgent task. Guided by the insight gained from two decades of modelization of 
spherically symmetric models, modern multidimensional calculations try to make a reli- 
able picture of the explosion on physical basis, minimizing the number of free parameters. 
Nevertheless, even supossing that the exploding object is a CO white dwarf near to the 
Chandrasekhar-mass limit the subject is not firmly resolved yet. In particular, there re- 
main uncertainties in the immediate pre-explosive stage and in the mechanism/s by which 
the thermonuclear flame is spreaded through the white dwarf. In this regard a difficult 
point is related with the acceleration of the combustion. We know that the explosion 
is triggered by a thermonuclear flame whose initial velocity is clearly subsonic. We also 
know, from parametrized calculations, that after several sound- crossing times the effective 
velocity of the deflagration must reach a significant fraction of of the local sound speed 
in order to synthesize the adequate amount of nuclear species to explain both light curve 
and spectra. This translates into an increase factor of thousands over the pure laminar 
value at maximum effective velocity. The physical mechanism behind this huge factor is 
currently attributed to flame surface corrugation by hydrodynamical instabilities (mainly 
Rayleigh- Taylor and turbulence). A branch in this scenario could be achieved if the al- 
ready fast deflagration turns into a steady detonation at late times, when flame density 
is below p ~ 4 10^ g.cm~^. This is the so called delayed-detonation scenario which, in 
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Table 1 

Main features of the moder 



t(s) 


Pc (g.cm 3) 


{pUame) (g-Cm ^) 




Mco{Mq) 


Ms,(Mo) 


MNi(M0) 





1.4 10^ 


1.06 10^ 


2.6 10-3 


1.29 


3.4 10"^ 


1.84 10^3 


1.54 


5.4 10^ 


2.04 10^ 


0.29 


0.65 


0.072 


0.272 


7.92 


3.5 10^ 




0.75 


0.34 


0.16 


0.54 



First column is the elapsed time since the beginning of the SPH simulation. Pre- 
viously, the initial stages of the explosion were calculated with a Lagrangian one- 
dimensional code. 



calculations carried out in spherical symmetry, have provided some of the best models for 
SNIa in the past |l]]. 

In this short paper we present the results of a calculation concerning the delayed deto- 
nation model carried out in three dimensions by using a SPH code. We describe the main 
features of the evolution of the model with more emphasis in the nucleoshynthetical side. 

2. Description of the model 

The initial model is a white dwarf of 1.36 Mq in hydrostatic equilibrium. The first stages 
of the explosion were followed by using a ID hydrocode until the central density declined 
to Pc = 4 10^ g.cm— 3. From here on the model was mapped into a 3D distribution of 
250,000 particules and its evolution followed with the SPH. In addition, the velocity field 
around the flame front was perturbed at the beginning in order to seed the hydrodynamical 
instabilities. This allows us to describe the large scale deformation of the flame in a self- 
consistent way although a lot of surface is lost owing to the limited resolution of the code. 
To take into account these hidden lengthscales a subgrid model has to be set. We assume 
that the main effect of the smallest lengthscales is to increase the laminar velocity of 
the flame, vi, until a first effective value fb > vi. Afterwards the final effective velocity 
Veff > t'b is directly computed from the hydrocode 0. 

During the progression of the calculation the geometrical features of the flame front were 
tracked by calculating the main scaling parameters of the surface such as Imin, Imax and 
flame fractal dimension D by using the method described in 0. As discussed in 0] the 
local evolution of D can be used as a practical criteria to turn the deflagration into a 
detonation. In particular, we have allowed such transition in those regions of the flame 
with fractal dimension higher than 2.5 once flame mean density is (puame) — 2 lO'' g.cm-^. 

The relevant features of the evolution of the model are summarized in Table 1. Basically 
the explosion proceeded in two phases. From t=0 s to t=1.54 s the combustion propagated 
subsonically as a deflagration. However, a considerable acceleration of the flame was 
seen because the effective flame velocity passed from v^fj — 8 lO'^c^ for t < 0.7s to 
Vcff — 0.14cs at t = 1.41s. At the end of this phase the total energy was already 
positive, although too low to represent a normal supernovae explosion. At t=1.54 s a 
detonation was artificially induced in those regions of the flame surface which had a large 
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Figure 1. Distribution of nuclear species Figure 2. Final production factors of scv- 
in velocity space eral isotopes normalized to solar system 

abundances 

fractal dimension. In consequence multiple detonations emerged from unconnected points 
spreaded on a very irregular surface. 

Information about the nucleosynthetic yields is provided in Table 1 and figures 1,2. 
The ejected amount of ^^Ni is enough to power the light curve of a typical Type la su- 
pernovae and the amount of silicon and other intermediate-mass elements is consistent 
with the spectrum near maximum light. A clear difference with respect ID calculations 
is that the abundance distribution in velocity space show a lot of dispersion in compari- 
son. The isotopic composition of the ejecta, calculated by postprocessing the dynamical 
model, is depicted in figure 2. Although we see a large production of vanadium there 
is not overproduction of problematic isotopes such as ^^Fe or ^^Ni. This is a feature 
of multidimensional models where electron captures proceed,on average, at lower density 
owing to the buoyancy of high entropy ashes. The ejected amount of carbon and oxygen is 
somewhat high, being distributed in isolated pockets within the star but not close to the 
center. Our results suggest that the three-dimensional version of the delayed detonation 
model is also able to reproduce the main observational features of Type la Supcrnovae. 

This work has been benefited from the MCYT grants EPS98-1348 and AYA2000-1785 
and by the DGES grant PB98-1183-C03-02. 
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